February 2, 2008 0:38 Proceedings Trim Size: 9in x 6in 



proc-SEWM 



PHASE CONVERSION AFTER A CHIRAL TRANSITION: 
EFFECTS FROM INHOMOGENEITIES AND FINITE SIZE' 



EDUARDO S. FRAG A 

Instituto de Fisica, Universidade Federal do Rio de Janeiro 
CP. 68528, Rio de Janeiro, RJ 21941-972, Brazil 
E-mail: fraga@if.ufrj.hr 



We discuss the kinetics of phase conversion, through the nucleation of bubbles and 
spinodal decomposition, after a chiral transition within an effective field theory 
approach to low-energy QCD. We study possible effects resulting from the finite 
size of the expanding system for both the initial and the late-stage growth of 
domains, as well as those effects due to inhomogeneities in the chiral field which 
act as a background for the fermionic motion. 



To model the mechanism of chiral symmetry breaking present in QCD, 
and to study the dynamics of phase conversion after a temperature-driven 
chiral transition, one can resort to low-energy effective models. In particu- 
lar, to study the mechanisms of bubble nucleation and spinodal decompo- 
sition in a hot expanding plasma, it is common to adopt the linear cr-model 
coupled to quarks, where the latter comprise the hydrodynamic degrees of 
freedom of the system^'^'^'*'^'^. The gas of quarks provides a thermal bath 
in which the long-wavelength modes of the chiral field evolve, and the lat- 
ter plays the role of an order parameter in a Landau-Ginzburg approach 
to the description of the chiral phase transition. The gas of quarks and 
anti-quarks is usually treated as a heat bath for the chiral field, with tem- 
perature T . The standard procedure is then integrating over the fermionic 
degrees of freedom, using a classical approximation for the chiral field, to 
obtain a formal expression for the thermodynamic potential of an infinite 
system. 

Let us consider a chiral field (j) — (a, tt ) , where cr is a scalar field and tt* 
are pseudoscalar fields playing the role of the pions, coupled to two flavors 
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of quarks according to the Lagrangian 

C = q[i^>^d^ + M,7° - M{4>)]q + \d^cj>d^'<t^ - V{ct^) . (1) 

Here q is the constituent-quark field q = {u, d) and fiq = ii/2> is the quark 
chemical potential. The interaction between the quarks and the chiral field 
is given by M((/)) = .g (cr + 275-? • tt) , and V{4i) = ^ (cr^ + tt^ - v^) - hga 
is the self-interaction potential for (f>. 

The parameters above are chosen such that chiral SUl{2) ® SUn{2) 
symmetry is spontaneously broken in the vacuum. The vacuum expectation 
values of the condensates arc (a) = f^c and (7?) = 0, where /^r = 93 McV is 
the pion decay constant. The explicit symmetry breaking term is due to 
the finite current-quark masses and is determined by the PCAC relation, 
giving hq = /tttoJ, where m.„ = 138 MeV is the pion mass. This yields 
v"^ = /|— mJ/A^. The value of = 20 leads to a cr-mass, to^ ~ 2A^/^+m^, 
equal to 600 MeV. For (/ > 0, the finite-temperature one-loop effective 
potential also includes a contribution from the quark fermionic determinant. 
In what follows, we treat the gas of quarks as a heat bath for the chiral 
field, with temperature T and baryon- chemical potential /i. Then, one 
can integrate over the fermionic degrees of freedom, obtaining an effective 
theory for the chiral field <f). Using a classical approximation for the chiral 
field, one obtains the thermodynamic potential 

fi(T,/z, = V{4>) - ^ IndetifG^i + W{4>)]IT} , (2) 

where Ge is the fermionic Euclidean propagator. From the thermodynamic 
potential one can obtain all the thermodynamic quantities of interest. 

To compute correlation functions and thermodynamic quantities, one 
has to evaluate the fermionic determinant within some approximation 
scheme. In ID systems one can usually resort to exact analytical methods'". 
In practice, however, the determinant is usually calculated to one-loop order 
assuming a homogeneous and static background chiral field. Nevertheless, 
for a system that is in the process of phase conversion after a chiral transi- 
tion, one expects inhomogeneities in the chiral field to play a role in driving 
the system to the true ground state. 

We propose an approximation procedure to evaluate the finite- 
temperature fermionic determinant in the presence of a chiral background 
field, which systematically incorporates effects from inhomogeneities in the 
chiral field through a derivative expansion. The method is valid for the case 
in which the chiral field varies smoothly, and allows one to extract infor- 
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mation from its long-wavclc^ngtli behavior, incorporating corrections order 
by order in the derivatives of the field. 

The Euler-Lagrange equation for static chiral field configurations 
contains a term which represents the fermionic density p{xo) = 
{vq/V) (xo + M{x))~^ \ S'o) , where |.to) is a position eigenstate with 

eigenvalue xq, and Uq = 12 is the color-spin- isospin degeneracy factor. In 
momentum representation: 



p{Xo)=VqTY, j ^ 



d?k_ 



-ik-x 



'y^iiojn + IJ') - 7 ■ k + M{x) 



(3) 



We can transfer the dependence to M{x) through a unitary transforma- 
tion, expand M{x + xo) around ^o, and use f ' = — iVfe. to write 



M-^ f d^k 1 

1 

1 + AM{-iVki,XQ) 



k + M{xo) 



(4) 



j°{i(j„ + IJ.) - J ■ k + M{xo) _ 

where xq is a c-number, not an operator, and AM{x,Xo) = WiM{xo)x^ + 
^ViVjM{xa)x'xi +■■■. 

If we focus on the long-wavelength properties of the chiral field and 
assume that the static background, M{x), varies smoothly and fermions 
transfer a small ammount of momentum to the chiral field, we can expand 
the expression above in a power series: 

d^k 



f d^k 



AM(-iVfe,,f) 



(27r)3 70(j^^ + ^) _ ^ . + M{x) 
1 



^°{iLOn + l^)-i-k + M{x) 



(5) 



AM{-iVk,,x) 



(6) 



which provides a systematic procedure to incorporate corrections brought 
about by inhomogeneities in the chiral field to the quark density, so that we 
can calculate p{x) = po{x) + pi{x) + P2ix) + ■ ■ ■ order by order in powers 
of the derivative of the background, M{x). The leading-order term in this 
gradient expansion for p{x) can be easily calculated and yields the well- 
known mean field result for p. The net effect of this leading term is to 
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correct the potential for the chiral field to Veff = V{^) + V^(0), where 



(7) 



where Ek{(j}) = \J k"^ + M{(f)). This sort of effective potential is commonly 
used as the thermodynamic potential in a phenomcnological description of 
the chiral transition for an expanding quark-gluon plasma created in a high- 
energy heavy-ion coUision^'^'^'^'^. However, the presence of a non-trivial 
background field configuration, e.g. a bubble, can in principle dramati- 
cally modify the Dirac spectrum, hence the determinant'''^. Results for the 
correction of the Laplacian term will be presented elsewhere^. 

In the process of phase conversion through bubble nucleation in a QGP 
of finite size, the set of all supercritical bubbles integrated over time will 
eventually drive the entire system to its true vacuum. The scales that de- 
termine the importance of finite-size effects are the typical linear size of 
the system, the radius of the critical bubble and the correlation length. 
For definiteness, let us assume our system is described by a coarse-grained 
Landau-Ginzburg potential, U (0, T), whose coefficients depend on the tem- 
perature. For the case to be considered, the scalar order parameter, ^, is 
not a conserved quantity, and its evolution is given by the time-dependent 
Landau-Ginzburg equation 

f =7[V^<^-f/UT)] , (8) 

where 7 is the response coefficient which defines a time scale for the system. 
The equation above is a standard reaction-diffusion equation, and describes 
the approach to equilibrium. 

If U {(j), T) is such that it allows for the existence of bubble solutions 
(taken to be spherical for simplicity), then supercritical (subcritical) bub- 
bles expand (shrink), in the thin- wall limit, with the following velocity: 

1 1 



(9) 



_Rc R{t) 

where Rc = {d— l)a/AF and AF is the difference in free energy between 

the two phases. The equation above relates the velocity of a domain wall 
to the local curvature. The response coefHcient, 7, can be related to some 
characteristic collision time. One can measure the importance of finite-size 
effects for the case of heavy-ion collisions by comparing, for instance, the 
asymptotic growth velocity (R » i?,c) for nucleated hadronic bubbles to 
the expansion velocity of the plasma. In the Bjorken picture, the typical 
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length scale of the expanding system is L(T) « {v^t,_>){Tc/T f = LoiTjTf, 
where is the collective fluid velocity and Lq = L(Tc) is the initial linear 
scale of the system for the nucleation process which starts at T < Tc- 

The relation between time and temperature provided by the cooling law 
that emerges from the Bjorken picture suggests the comparison between the 
following "velocities" : 

..-^^-f^USVVi-^l . (10) 



dT \2v^aT^ J \ T J \ T, , 

the asymptotic bubble growth "velocity", and the plasma expansion "ve- 
locity" vl = (dL/dT) = -{3Lo/Tc){Tc/T)'^. The quantity 6 is a number of 
order one to first approximation, and comes about in the estimate of the 
phenomenological response coefficient 7(T) « b/2T. Using the numerical 
values adopted previously and a/T^' ~ 0.1, we obtain-^" 

"'-^m-il . (11) 



VL v,\T 

One thus observes that the bubble growth velocity becomes larger than 
the expansion velocity for a supercooling of order « Vz/20 < 5%. A 
simple estimate points to a critical radius larger than 1 fm at such values 
of supercooling**. Therefore, finite-size effects appear to be an important 
ingredient in the phase conversion process right from the start in the case 
of high-energy heavy-ion collisions'*^. 
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